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Abstract 

We investigate the spreading of passive tracers in closed basins. If the charac¬ 
teristic length scale of the Eulerian velocities is not very small compared with 
the size of the basin the usual diffusion coefficient does not give any relevant 
information about the mechanism of spreading. 

We introduce a finite size characteristic time t{6) which describes the 
diffusive process at scale 5. When 5 is small compared with the typical length 
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of the velocity field one has t{5) ~ A“^, where A is the maximum Lyapunov 
exponent of the Lagrangian motion. At large 5 the behavior of t{S) depends 
on the details of the system, in particular the presence of boundaries, and in 
this limit we have found a universal behavior for a large class of system under 
rather general hypothesis. 

The method of working at fixed scale S makes more physical sense than the 
traditional way of looking at the relative diffusion at fixed delay times. This 
technique is displayed in a series of numerical experiments in simple ffows. 
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I. INTRODUCTION 


The understanding of diffusion and transport of passive tracers in a given velocity field 
has both theoretical and practical relevance in many fields of science and engineering, e.g. 
mass and heat transport in geophysical flows (for a review see |]^,|^), combustion and chemical 
engineering [^. 

One common interest is the study of the mechanisms which lead to transport enhance¬ 
ment as a fluid is driven farther from the motionless state. This is related to the fact that 
the Lagrangian motion of individual tracers can be rather complex even in simple laminar 


flows iJH- 


The dispersion of passive scalars in a given velocity field is the result, usually highly non¬ 
trivial, of two different contributions: molecular diffusion and advection. In particular, one 
can have rather fast transport, even without molecular diffusion, in presence of Lagrangian 
chaos , which is the sensitivity to initial conditions of Lagrangian trajectories. In addition, 
also for a 2D stationary velocity field, where one cannot have Lagrangian chaos 0 , in pres¬ 
ence of a particular geometry of the streamlines the diffusion can be much larger than the 
one due only to the molecular contribution, as in the case of spatially periodic stationary 


flows 0J|]. 


Taking into account the molecular diffusion, the motion of a test particle (the tracer) is 
described by the following Langevin equation: 



( 1 ) 


where u(x, t) is the Eulerian incompressible velocity field at the point x and time t, r]{t) is 
a Gaussian white noise with zero mean and 


< Vi{t)Vj{t) >= ‘2L)ohijS{t - t ), 


( 2 ) 


where Dq is the (bare) molecular diffusivity. 

Denoting 0(x, t) the concentration of tracers, one has: 
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dtQ + {u-V)Q = DqAQ. 


( 3 ) 


For an Eulerian velocity field periodic in space, or anyway defined in infinite domains, 
the long-time, large-distance behavior of the diffnsion process is described by the effective 
diffnsion tensor Df- {eddy-dijfusivity tensor): 


Dfj = lim — < (a;i(t)- < Xi >){xj{t)- < xj >) > , 
■' t^oo Zt 


( 4 ) 


where now x(t) is the position of the the tracer at time t, i,j = 1, - ■ ■ ,d (being d the spatial 
dimension) , and the average is taken over the initial positions or, eqnivalently, over an 
ensemble of test particles. The tensor Dfj gives the long-time, large-distance eqnation for 
< 0 > i.e. the concentration field locally averaged over a volnme of linear distance mnch 
larger than the typical length of the velocity field, according to 


d q2 


< 0 > . 


( 5 ) 


The above case, with hnite is the typical sitnation where the diffnsion, for very large 
times, is a standard diffnsion process. However there are also cases showing the so-called 
anomalous diffusion: the spreading of the particles does not behave linearly with time bnt 
has a power law t‘^'^ with u ^ 1/2. Transport anomalies are, in general, indicators of the 
presence of strong correlation in the dynamics, even at large time and space scales 

In the case of inhnite spatial domains and periodic Enlerian helds the powerfnl mnltiscale 
techniqne (also known as homogenization in mathematical literatnre) gives a nsefnl tool for 
stndying standard diffnsion and, with some precantions, also the anomalons sitnations . 

On the other hand we have to stress the fact that diffnsivity tensor (|^) is mathematically 
well defined only in the limit of inhnite times, therefore it gives a sensible resnlt only if the 
characteristic length lu of the velocity held is mnch smaller than the size L of the domain. 

The case when and L are not well separated is rather common in many geophysical 
problems, e.g. spreading of pollntants in Mediterranean or Baltic sea, and also in plasma 
physics. Therefore it is important to introdnce some other characterizations of the dihnsion 
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properties which can be used also in non ideal cases. For instance, propose to employ 
exit times for the study of transport in basins with complicated geometry. 

In Section we introduce a characterization of the diffusion behavior in terms of the 
typical time t{6) at scale 5; this allows us to define a finite size diffusion coefficient D{6) ~ 
/^From the shape of t{6) as a function of 6, one can distinguish different spreading 

regimes. 


In Section |I| we present the results of numerical experiments in closed basins and present 
new results relative to the behavior of the diffusion coefficient near the boundary (a detailed 
discussion is in the appendix). 

In Section 0 we summarize our results and present conclusions and we discuss the 
possibility of treatment of experimental data according to the method introduced in Section 

0 


II. FINITE SIZE DIFFUSION COEFFICIENT 


Before a general discussion let us start with a simple example. Consider the relative 
diffusion of a cloud of N test particles in a smooth, spatially periodic velocity field with 
characteristic length /„. We assume that the Lagrangian motion is chaotic i.e. the maximum 
Lyapunov exponent A is positive. Denoting with R‘^(t) the square of the typical radius of 
the cloud 


R^{t) =-C |xi(t)- -C Xi(t) > 1^ > , 


where 


1 ^ 

< x*(f) >= 

i=i 


we expect the following regimes to hold 


-- 2 ^ ^^(0)exp(L(2)t) 

2Dt 


if 


( 6 ) 

(7) 

( 8 ) 


5 









where L{2) > 2A is the generalized Lyapunov exponent |1l2|,p!3|, D is the diffusion coefficient 
and the overbar denotes the average over initial conditions. 

In this paper we prefer to study the relative diffusion instead of the usual absolute 
diffusion. For spatially infinite cases, without mean drift there is no difference; for closed 
basins the relative dispersion is, for many aspects, more interesting than the absolute one 
and, in addition, the latter is dominated by the sweeping induced by large scale flow. 

Furthermore we underline that although the dynamics of the ocean circulation is domi¬ 
nated by large mesoscale gyres, the smaller scales activities within the gyres control impor¬ 
tant local phenomena as deep water formation in North Atlantic and in Mediterranean basin 


141. Therefore the study of relative diffusion could be relevant to describe this small-scale 


motion and can give crucial informations on the way to parameterize the subgrid scales in 


ocean numerical global model |15 


Another, at first sight rather artificial, way to describe the above behavior is by in¬ 
troducing the “doubling time” r(5) at scale 6 as follows: we define a series of thresholds 
^(n) _ ^ where is the initial size of the cloud, defined according to (^, and then 

we measure the time it takes for the growth from to 5*-^^ = rS^^\ and so on for 

r(5d)), 

up to the largest scale under consideration. For the threshold rate r 
any value can be chosen but too large ones might not separate different scale contributions, 
though strictly speaking the term “doubling time” refers to the threshold rate r = 2. 

Performing A/” 3> 1 experiments with different initial conditions for the cloud, we define 
the typical doubling time t{S) at scale 6 as 


rW =< T(S) 

i=l 


(9) 


Let us stress the fact that the average in is different from the usual time average. 

From the average doubling time we can define the finite size Lagrangian Lyapunov ex¬ 
ponent as 

Inr 


A(5) = 


r(5)’ 


( 10 ) 


6 




















which is a measure of the average rate of separation of two particles at a distance 6. Let us 
remark that A(5) is independent of r, for r —> !■*■. For very small separations (i.e. 6 <C lu) 
one recovers the standard Lagrangian Lyapunov exponent A, 

A = lim—— Inr. (11) 

<5^0 r(5) ^ ' 

See for a detailed discussion about these points. In this framework the hnite size diffusion 
coefficient D{6) dimensionally turns out to be 

D{S) = S^X{6). (12) 


Note the absence of the factor 2, as one can expect by the dehnition (^, in the denominator 
of D{6) in equation (p)2D; this is due to the fact that t(6) is a difference of times. For a 
standard diffusion process -D(S) approaches the diffusion coefficient D (see eq. (||)) in the 
limit of very large separations (6 3> lu)- This result stems from the scaling of the doubling 
times t(6) ~ for normal diffusion. 

Thus the hnite size Lagrangian Lyapunov exponent A(5), or its counterpart D(S), embody 
the asymptotic behaviors 


A(i) 


(13) 


A if S lu 
if 5 > 

One could naively conclude, matching the behaviors at 5 ~ lu, that D ~ A/^. This is not 
always true, since one can have a rather large range for the crossover due to the fact that 
nontrivial correlations can be present in the Lagrangian dynamics 0- 


Another case where the behavior of t{6) as a function of 6 is essentially well understood 
is 3D fully developed turbulence. For sake of simplicity we neglect intermittency effects. 
There are then three different ranges: 


1. 6 rj = Kolmogorov length : l/r((5) ~ A; 

2. 7] -C 5 -C / = typical size of the energy containing eddies: from the Richardson law 
R2(t) ~ one has l/r((5) ~ 
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3. 5 S> / : usual diffusion behavior l/r(5) ~ <5 ^ . 


One might wonder that the proposal to introduce the time t{6) is just another way to 


look at as a function of t. This is true only in limiting cases, when the different 

characteristic lengths are well separated and intermittency is weak. In |]T8|-pm rather close 
techniques are used for the computation of the diffusion coefficient in nontrivial cases. 

The method of working at hxed scale 6, allows us to extract the physical information at 
that spatial scale avoiding unpleasant troubles of the method of working at a hxed delay 
time t. For instance, if one has a strong intermittency, and this is a rather usual situation, 
R^{t) as a function of t can appear very different in each realization. Typically one can 
have, see hgure ^a, different exponential rates of growth for different realizations, producing 


a rather odd behavior of the average R‘^{t) without any physical meaning. For instance in 


hgure mb we show the average R^{t) versus time t; at large times we recover the dihusive 
behavior but at intermediate times appears an apparent anomalous regime which is only due 
to the superposition of exponential and dihusive contributions by diherent samples at the 
same time. On the other hand exploiting the tool of doubling times one has an unambiguous 
result (see hgure ^c). 

Of course the interesting situations are those where the diherent characteristic lengths 


{rj, /, L) are not very diherent and therefore each scaling regime for R‘^{t) is not well 
evident. 


III. NUMERICAL RESULTS 


Here we present some numerical experiments in simple models with Lagrangian chaos 
in the zero molecular dihusion limit. Before showing the results, we describe the numerical 
method adopted. 

We choose a passive tracers trajectory having a chaotic behavior, i.e. with a positive 
maximum Lyapunov exponent, computed by using standard algorithms |^|. Then we place 
N — 1 passive tracers around the hrst one in a cloud of initial size 
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R{0) = 5(0) = 5 ( 0 ), 


with -R(O) defined by equation (j^). In order to have average properties we repeat this 
procedure reconstructing the passive cloud around the last position reached by the reference 
chaotic tracer in the previous expansion. This ensures that the initial expansion of the cloud 
is exponential in time, with typical exponential rate equal to the Lyapunov exponent. 

Further we define a series of thresholds (as described in Section 2) n = 

'n^max and we measure the time T„ spent in expanding from to . The value 

of rijnax has to be chosen in such a way that ~ ^max-, where 5max corresponds to 

the uniform distribution of the tracers in the basin (see forthcoming discussion and the 
Appendix). Each realization stops when 5{t) = 

Therefore following |jT^ we define a scale dependent Lagrangian Lyapunov exponent as: 




Inr = 


In r. 


(14) 


<Tn>e'~~ r(5W) 

In equation (p!4D we have implicitly assumed that the evolution of the size 6{t) of the cloud 
is continuous in time. This is not true in the case of discontinuous processes such as maps 
or in the analysis of experimental data taken at fixed delay times. Denoting T„ the time 
to reach size S > from , now 5 is a fluctuating quantity, equation (JT^) has to be 


modified as follows 16 


A(5(’")) = 




In 


5W 


( 16 ) 


In our numerical experiments we have the regimes described in sect. 2: exponential 
regime , i.e. A(5) = A, and diffusion-like regime i.e. A(5) = D/5^, at least if the size L of 
the basin is large enough. 

For cloud sizes close to the saturation value dmax we expect the following behavior to 
hold for a broad class of systems: 


A(5) 


n(d) (Smax-S) 

5 


(16) 
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The constant of proportionality is given by the second eigenvalne of the Perron-Frobenins 
operator which is related to the typical time of exponential relaxation of tracers’ density 
to the nniform distribntion Actnally, the analytical evalnation of this eigenvalne can be 
performed only for extremely simple dynamical systems (for instance random walkers, as 
shown in the Appendix). As a conseqnence the range of validity for (P!B|) can be assessed 
only by nnmerical simnlation. 


A. A model for transport in Rayleigh-Benard convection 


The advection in two dimensional incompressible flows is described, in absence of molec- 
nlar diffnsion, by Hamiltonian eqnation of motion where the Hamilton fnnction is the stream 
fnnction ip: 


dx dip dy dip 

dt dy ’ dt dx 


(17) 


If Ip is time-dependent the system ([T7|) is non-antonomons and in general non-integrable, 
then chaotic trajectories may exist. 

One example is the model introdnced in to describe the chaotic advection in the 
time-periodic Rayleigh-Benard convection. It is dehned by the stream fnnction: 


ip{x, y,t) = — sm{k[x + B sin(a;t)] }W{y), 

K 


(18) 


where Wpy) is a fnnction that satishes rigid bonndary conditions on the snrfaces y = 0 and 
y = a (we nse W{y) = sm{7iy/a)). The direction y is identihed with the vertical direction 
and the two snrfaces y = a and y = 0 are the top and bottom snrfaces of the convection cell. 
The time dependent term B sin(a;t) represents lateral oscillations of the roll pattern which 
mimic the even oscillatory instability ||22|| . 

Trajectories starting near the roll separatrices conld have positive Lyapnnov exponent 
and thns display chaotic motion and diffnsion in the x direction. It is remarkable that in 
spite of the simplicity of the model, the agreement of the nnmerical resnlts with experimental 
ones is qnite good 0. 
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Defining a passive cloud in the x direction (i.e. a segment) and performing the expansion 
experiment described in the previous section we have that, until 5 is below a fraction of the 
dimension of the cell, A(5) = A (hgure |^a). For larger values of 5 we have the standard 
diffusion A((5) = with good quantitative agreement with the value of the diffusion 

coefficient evaluated by the standard technique, i.e. using E?{t) as a function of time t 
(compare hgure ^a with hgure Hb). 

To conhne the motion of tracers in a closed domain, i.e. x G [—L,L], we must slightly 
modify the streamfunction (|^. We have modulated the oscillating term in such a way that 
for |a;| = L the amplitude of the oscillation is zero, i.e. B —>■ Bsin{nx/L) with L = 2'xn/k 
{n is the number of convective cells). In this way the motion is conhned in [—T, L\. 

In hgure ^ we show X{S) for two values of L. If L is large enough one can well see the 
three regimes, the exponential one, the dihusive one and the saturation given by equation 
(pT|). Decreasing L the range of the dihusive regime decreases, and for small values of L it 
disappears. 


B. Modified Standard Map 


One of the simplest deterministic dynamical system displaying both exponential growth 
of separation for close trajectories and asymptotic dihusive behavior is the standard 
(Chirikov - Taylor) mapping |Q. It is customarily dehned as 

Xn+1= Xn + Ksinyn 
Vn+l = yn + Xn+1 mod 271 

This mapping conserves the area in the phase space. It is widely known that for large enough 
values of the nonlinearity strength parameter iF S> iFc — 1 the motion is strongly chaotic 
in almost all the phase space. In this case the standard map, in the x-direction mimics the 
behavior of a one-dimensional random walker, still being deterministic, and so one expects 
the behavior of A(5) to be quite similar to the one already encountered in the model for 


Rayleigh-Benard convection without boundaries. Numerical iteration of (jT^) for a cloud of 
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particles clearly shows the two regimes described in ([T^), similar to that showed for the 
model discussed in the previous section. 


We turn now to the more interesting case in which the domain is limited by boundaries 
reflecting back the particle. To achieve the conhnement of the trajectory inside a bounded 
region we modify the standard map in the following way 


{ Xn+l =Xn + Kf{Xn+l) sini/n 

Vn+l = yn + Xn+l “ Kf'{Xn+l) COS mod 271. 


( 20 ) 


where f{x) is a function which has its only zeros in ±L. Since the mapping is dehned in 
implicit form, the shape of / must be chosen in such a way to assure a unique dehnition for 


{xn+i,yn+i) given {xn,yn)- For any / fulhlling this request the mapping (poD conserves the 
area. A trial choice could be 


fix) = 


1 

L 


\x\ 


|x| < t 

I <\x\ < L 


( 21 ) 


L-e 

Strictly speaking this is not quite an appropriate choice, since it renders the map discon¬ 
tinuous at |x| = i, but this is an irrelevant point and it is easy to bypass this obstacle by 


assuming a suitably smoothed version of (pT]). 

Performing the doubling times computation (P) one recovers both the exponential and 
diffusive regimes for A (5), and in addition one has the saturation regime (^). Figure ^ 
shows the behavior of the scale dependent diffusion coefficient D{6) (|T^. Approaching the 
saturation value 6max the diffusion coefficient quickly drops to zero, following the asymptotic 


law ([Tq) derived in the appendix. The qualitative behaviors in hgure ^ do not depend on 
the details of the function /. 


C. Point vortices in a Disk 


As another example, we consider the two-dimensional time-dependent flow generated by 


the motion of N point vortices in a closed domain ||^. For a unitary disk the positions 
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of the vortices {xi = r* cos 9i,yi = rjsin 6*j), with circulation Fj, evolve according to the 
Hamiltonian dynamics 

I dH . _ I dH 

Fj dyt ’ Fj dxi 

where the Hamiltonian is 



" = -i^i:nniog 

i>j 


r‘f + r‘j — cos(9i — 9j) 
1 + r?r| 


N 


2rirj cos{9i — 9j 
Passive tracers evolve according to m with "0 given by 

+ rf — 2rri cos(6* — 6*j) 




(23) 


1 " 

iix,y) = -T-XJmog 


1 + r'^rf — 2rri cos{9 — 9i) 


(24) 


where {x = r cos 9,y = r sin 9) denote the tracer position. 

Figure ^ shows the relative diffusion as a function of time in a system with 4 vortices. 
Apparently there is an intermediate regime of anomalous diffusion. On the other hand from 
hgure 1^ one can see rather clearly that, with the method of working at hxed scale, only two 
regimes survive: the exponential one and that one due to the saturation. Comparing hgure 
P and hgure ^ one understands that the mechanism described in Section 2 has to be held for 
responsible of this spurious anomalous dihusion. We stress the fact that these misleading 
behaviors are due to the superposition of diherent regimes and that the method of working 
at hxed scale has the advantage to eliminate this trouble. 

The absence of the dihusive range A(^) ~ 5“^ is due to the fact that the characteristic 
length of the velocity held, which is comparable with the typical distance between two close 
vortices, is not much smaller than the size of the basin. 


IV. CONCLUSIONS 

In this paper we investigated the relative dispersion of passive tracers in closed basins. 
Instead of the customary approach based on the average size of the cloud of tracers as a 
function of time, we introduced a typical inverse time A (5) which characterizes the dihusive 
process at hxed scale 6. 
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For very small values of S, A(5) coincides with the maximum Lagrangian Lyapunov 
exponent which is positive in the case of chaotic Lagrangian motion. For larger 6 the shape 
of A(5) depends on the detailed mechanism of spreading which is given by the strncture of 
the advecting flow, which is in turn conditioned by the presence of boundaries. In the case of 
diffnsive regime, one expects the scaling A((5) ~ S~‘^, which leads to a natural generalization 
of the diffnsion coefficient as D{6) = X{6)6‘^. 

The effectiveness of finite size quantities A(6) or D(S) in characterizing the dispersion 
properties of a cloud of particles is demonstrated by several numerical examples. 

Furthermore, when 5 gets close to its satnration valne (i.e. the characteristic size of the 
basin), a simple argument gives the shape of A(5) which is expected to be universal with 
respect to a wide class of dynamical systems. 

In the limiting case when the characteristic length of the Eulerian velocity lu and the 
size of the basin L are well separated, the customary approach and the proposed method 
give the same information. In presence of strongly intermittent Lagrangian motion, or when 
lu/L is not much smaller than one, the traditional method can give misleading results, for 
instance apparent anomalous scaling over a rather wide time interval, as demonstrated by 
a simple example. 

We want to stress out that our method is very powerful in separating the different scales 
acting on diffusion and consequently it could give improvement abont the parameterization 
of small-scale motions of complex flows. The proposed method could be also relevant in the 
analysis of drifter experimental data or in numerical models for Lagrangian transport, in 
particular for addressing the question abont the existence of low dimensional chaotic flows. 
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APPENDIX 


In this appendix we present the derivation of the asymptotic behavior (P^bf) of A(5) for 5 
near to the saturation, for a one dimensional Brownian motion in the domain [—P, P], with 
reflecting boundary conditions. The evolution of the probability density p is ruled by the 
Fokker-Planck equation 


dp 1 ^ d'^p 
dt 2 dx"^ 


(25) 


with the Neumann boundary conditions 

|(±L) = 0, (26) 

The general solution of (^) is 

CXD 

p{x,t)= ^ + C.C (27) 

k=—oo 

where 

, fc = 0 ,±l,±2,... (28) 

At large times p approaches the uniform solution po = 1/2L. Writing p as p{x,t) = Po + 
6p{x, t) we have, for f S> ri , 


6p ~ exp(—f/ri). 

The asymptotic behavior for the relative dispersion R‘^{t) is 

R^it) ^ ^ y ~ x'Yp{x^ t) dx dx' 

For f 3> Ti using (P9[) we obtain 


(29) 


( 30 ) 
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^ i Ae-'/’' 


( 31 ) 


Therefore for 5{t) = R{t) one has 




(32) 


The saturation value of 6 is 6max = -^/\/3; so for f S> ri, or equivalently for {6max — S)/S I, 
we expect 


d , r \ /■ r\ 1 ^max ^ 

— Ind = Aid) =-;- 

dt ^ ^ Ti 6 


(33) 


which is (pTl). 

Let us remark that in the previous argument for A (5) for S ~ Smax the crucial point is the 
exponential relaxation to the asymptotic uniform distribution. In a generic deterministic 
chaotic system it is not possible to prove this property in a rigorous way. Nevertheless one 
can expect that this request is fulhlled at least in non-pathological cases. In the terminology 
of chaotic systems the exponential relaxation to asymptotic distribution corresponds to have 
the second eigenvalue a of the Perron-Frobenius operator inside the unitary circle; now the 


relaxation time is ti = — In I a I p5[ . 
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FIGURES 

FIG. 1. a) Three realizations of as a function of t built as follows; R‘^{t) = exp( 27 f) 


if R^{t) < 1 and R^{t) = 2D{t — f*) with 7 = 0.08,0.05,0.3 and do = 10“^, D = 1.5. b) R?‘{t) 
as function of t averaged on the three realizations shown in figure la. The apparent anomalous 
regime and the diffusive one are shown, c) A(d) vs 5, with Lyapunov and diffusive regimes. 

FIG. 2. Lagrangian motion given by the Rayleigh-Benard convection model with: 
A = 0.2, B = 0.4, uj = 0.4, k = 1.0, a = vr, the number of realizations is Af = 2000 and 
the series of thresholds is 5n = dor"" with do = 10“^ and r = 1.05. a)A(d) vs d, the horizontal line 
indicates the Lyapunov exponent A = 0.022, the dashed line is D6~^ with D = 0.26. h)R'^{t) as 
function of t, the line is 2Dt with D = 0.26. 

FIG. 3. A(d) vs d for the same model and parameters of figure 2 , but in a closed domain with 
6 (crosses) and 12 (diamonds) convective cells. The lines are respectively: (a) Lyapunov regime 
with A = 0.017; (b) diffusive regime with D = 0.021; (c) saturation regime with 5max = 19.7; (d) 
saturation regime with 5max = 5.7. 

FIG. 4. D{5) vs d for the modified standard map with K = 8, L = 1000 and I = 990. The 

series of thresholds is 5n = dor"" with do = 10“^ and r = 2^/^®. The horizontal line indicates 
the diffusion coefficient in the limit of infinite system, the dashed curve represents the saturation 
regime. 

FIG. 5. R‘^{t) for the four vortex system with Fi = F 2 = —Fa = —r 4 = 1. The threshold 

parameter is r = 1.03 and do = 10“^, the dashed line is the power law R‘^{t) ~ The number 
of realizations is A7 = 2000 . 
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FIG. 6. A((5) vs 5 for the same model and parameters of figure 5. The horizontal line indicates 


the Lyapunov exponent (A = 0.14), the dashed curve is the saturation regime with 5max = 0.76. 
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